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. ABSTRACT 

Previous work on parameter plane, three dimensional pa- 
rameter space, and singular lines in the parameter plane are 
reviewed. 

A general concept of n-dimensional parameter space is 
hypothesized whereby the parameter plane becomes a special 
case of the general hypothesis. By the same argument the 
singular line is shown to be a special case of the singular 
hyperplane . 

Existence criteria for singular lines are established, 
and compensation methods for creating singular lines in non- 
singular systems are derived and used. 
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CHAPTER I 



INTRODUCTION 



The analysis of linear control systems by use of the 
Laplace Transforms of the systems' characteristic differential 
equations is adequately accomplished by root locus or fre- 
quency response methods when all coefficients of the complex 
variable s are real constants. 

These same methods are also sufficient to effect anal- 
ysis when the coefficients are functions of a single real var- 
iable K. The general characteristic equation 

n , 

f ( s ) = l a. s K = 0 (1-1) 

k=0 K 

may be expressed as 

n , 

f ( s ) = l [b, (K)+c. )s K = 0 , (1-2) 

k=0 k K 

from which follows 



I 

k=0 

n 



I 

k=0 



b k (K)s J 



c, s 
k 



-1 , 



(1-3) 



which is then in the proper form for either root locus or fre- 
quency response analysis. 

However, many modern control systems, although still lin- 
early conceived or approximated, cannot be handled by the 
above somewhat "classical" techniques. The complexities 
accompanying the advancing "state of the art" in control system 
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and compensation design often result in characteristic equation 
coefficients that are functions of two or more variables. 

In 1959, Mitrovic introduced a method [1] for determining 
the roots of the characteristic equation when the coefficients 
of the two lowest ordered terms vary independently. 

A point in the s-plane, 

s = -a + jco (1-4) 

may be specified in polar form as 

s = -£co + jco /i-C 2 , (1-5) 

where co n is the radial distance from the origin to the point 
and £ is the cosine of 0, the angle between the real axis and 

0) . 

n 

Since 



£co = COS0 
n 



and 



co /l-C 2 = sin0, 
n ' 



equation 1-5 may be expressed as 



( 1 - 6 ) 



(1-7) 



and 



s = to n (cos0+jsin0) ; 



s^ = co k (cosk0 + jsink0) 
n 



(1-8) 



(1-9) 



As previously defined, for left-hand plane roots 



0 



-1 

cos 



(-£) ; 



( 1 - 10 ) 



therefore 



cosk0 = cosk[cos ^(-?)], (1-11) 
which can be expressed as a Chebyshev function of the first 
kind , T k (-£) , and 

sink© = sink[cos ^(-£)], (1-12) 
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+Im 




Fig. 1.1 A Point in the s-Plane 
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which can be expressed by a Chebyshev function of the second 



kind, s in 0 U k ( - ? ) . This permits equation 1-9 to be written 
s k = aj k [T k (-«;)+jsin0U k (-^) ] (1-13) 

s k = u> k [(-l) k T k U)+jsin0 (-l) k+1 U k U) ] (1-14) 

and equation 1-1 to be written 

f ( s ) = l a, coN-l) T, (?)+jsin0 £ a,/(-l) U,(?) = °« 

k=0 K n K k=0 K n K 

(1-15) 



Since, for any non-zero value of 0, both the real and imag- 
inary parts of equation 1-15 must independently equal zero, 
two simultaneous equations for a k are formed: 

k , , , k„ 



n 



l a ojN- 1) k T (?) = 0 
k=0 * n K 



(1-16) 



l a. a) k (-l) k+1 U, (?) = 0 . (1-17) 

k=0 K n K 

The following relationship exists between the two kinds of 
Chebyshev functions : 

T k U) = CU k (C) " U k-l U) ; d-3-8) 

also 

£ k+1 k 5 V k 

I (-1) a A (?) = - l (-1)\0.(0; d-19) 

k=0 k n k k=Q k n k 

therefore for any non-zero value of ?, the simultaneous equa- 
tions reduce to 

X '-"VX- 1 ( '> - 0 (1 - 20) 

k=0 



l (-l) k a, 03 k U, (?) = 0 . (1-21) 

k=0 k n K 
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U k (£) can be obtained for all values of k from the recur- 



sion formula 

U k+l (c) = 2cU k (C) ” U k-l (c)/ 
where, for all values of £ , 

U_! ( C) = -1 
U Q (C) = 0 

U^C) = 1 • 



Equation 1-1 can be rewritten as 

f(s) = a s n + a , s n ^ + . . . a~s 2 + A, s + A n 
n n-1 2 1 U 

and equations 1-20 and 1-21 become 



-A, 



n 

+ l 

k=2 



(-l) k a, w k U, , U) = 
k n k-1 



( 1 - 22 ) 

(1-23) 

(1-24) 

(1-25) 

0, (1-26) 
(1-27) 



-“n A l + J, " 0 ' (1 - 28) 

k=2 

Thus, Aq and become the coordinate axes of the 
Mitrovic Plane; and lines of constant £ and o> n , along with 
constant real root lines from the direct substitution of (-a) 
for s in equation 1-26, can be plotted, giving the solution 
of the characteristic equation for any point (Aq , A^) . 

This method was later extended by Elliott, Thaler, and 
Heseltine [2,3] for specific coefficient pairs, and then gen- 
eralized by Siljak [2,4] for any pair of variable coefficients 
in the Coefficient Plane method, where equation 1-1, when 
written as 

f(s) = A, s 1 + A s m + T a,s k = 0, (1-29) 

1 m k=0 K 

k^l ,m 
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results in the set of simultaneous equations in (A.^, A^) 





k^l ,m 



(1-30) 





k=0 




0. 



k^l ,m 



(1-31) 



Although this is an improvement of Mitrovic's original 
method, it is still limited by the requirement that only two 
coefficients of the characteristic equation contain variable 
elements. However, in 1964 Siljak [5] extended this process 
into the Parameter Plane, where two variable parameters can 
appear in any of the coefficients of the characteristic equa- 
tion. In the linear case. 



where a and 3 are the adjustable parameters. Hence, equation 
1-1 becomes 



a k ■ b k° + c k 6 + f k- 



(1-32) 



K=U K=U K=U 

and equations 1-20 and 1-21 become 




n 




(?) + 6 l (-l) k c 



k=0 





k=0 




k=0 k=0 




k=0 
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For notational ease. 




k=0 



(1-36) 



(1-37) 



and 



B^ + + F 1 = 0 



(1-38) 



B 2 a + C 2 B + F 2 = 0 



(1-39) 



which are solved by standard methods, giving constant £, to n , 
and a lines in the a-0 Parameter Plane for the solution of the 
characteristic equation. 

In 1968 Cadena [6] extended the parameter plane concept 
to three independent parameters, where 



which when solved give the line of intersection of two planes 
in 3-parameter space. By an ingenious device of rotating the 
y coordinate by 90° to form the (-3) ordinate in the a-0 plane, 
Cadena, by graphical projections, reduces the constant t; , io , 
and 0 planes to families of lines in the' a-0 plane — each 
line representing a specific value of y. In the same work, it 
is shown that for an n-parameter system, n-roots may be fixed 
by the adjustment of each parameter to a specific value. In 
this case the remaining roots are then determined by factoring 



a k = b k a + c k S + d k Y + f k . 



(1-40) 



The parametric simultaneous equations then become 



B^a + C^0 + D^y + F^ = 0 



(1-41) 



B 2 a + C 2 0 + D 2 y + F 2 = 0 



(1-42) 
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the quotient which results from the division of the character- 
istic equation by the product of the specified roots. 

What emerges from this "reduced characteristic equation" 
concept is that a designer has at his disposal a "degree of 
freedom" for each variable parameter in his system. If, for 
example, the system contains three parameters , three roots 
may be fixed by rigid adjustment of all three parameters; if 
only two roots are specified (e.g., a complex conjugate pair), 
the reduced characteristic equation then becomes a function of 
one parameter (the remaining degree of freedom) , and ordinary 
root locus methods determine the remaining roots; if only one 
real root is specified, the reduced characteristic equation 
is then a two parameter case , leading to root solution by 
Parameter Plane methods. 

Thus Parameter Plane techniques have become a powerful 
tool for both the analysis and synthesis of modern complex 
control systems. 
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CHAPTER II 



THE SINGULAR LINE 



In the preceding chapter it was shown that constant Z, , u) n , 
and a lines could be mapped in the a-B plane by the simulta- 
neous solutions of equations 1-38 and 1-39. Using Cramer's 
Rule, points M^(£,u) n ) can be mapped into a point M(a,B) by 



and 



a = 



' F 1 C 1 



' F 2 C 2 
B 1 C 1 



B 2 C 2 



C 1 F 2 ~ C 2 F 1 
B 1 C 2 “ B 2 C 1 



6 = 



B 1 - F 1 



B 2 ~ F 2 
B 1 C 1 



B 2 C 2 



B 2 F 1 



B 1 F 2 



B 1 C 2 



B 2 C 1 



( 2 - 1 ) 



( 2 - 2 ) 



the point, M, being the intersection of constant £ lines and 
constant u> lines. By this method the roots of characteristic 
equations can be determined from the output graphs of various 
Parameter Plane computer programs, using linear interpolation, 
if necessary, to find all roots at any point M(a,8). 

An obvious limitation of this method is that the deter- 
minant of the coefficient matrix be non-singular, or non-zero. 

Although this limitation existed, it was dismissed as be- 
ing of little significance until 1967, when Bowie [7], in 
attempting to solve a sixth order characteristic equation. 
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found, that by selecting various points M(a,8), only certain 
roots could be found using Parameter Plane methods. However, 
by substituting the co-ordinates of the point M(a,6) into the 
characteristic equation and factoring, the undetermined roots 
were found to be complex conjugate pairs. Moreover, many 
points M(a,$) were found to have the same complex pair in com- 
mon; these common root pair points formed a straight line in 
the a-3 plane which was the locus of a constant C-constant o> n 
pair. 

Upon further investigation of this phenomenon, Bowie dis- 
covered that lines of constant values of £, to n , which he called 
singular lines, existed whenever equations 1-38 and 1-39 were 
linearly dependent. 

If the coefficient matrix of equations 1-38 and 1-39 is 
expanded in terms of £ and a) , an infinite number of singular 
lines can be found by setting the determinant to zero and solv- 
ing for either £ or w in terms of the other. For an equation 
of order n, there are a maximum of 2(n-l) values of which 
produce singular lines for each of the infinite number of val- 
ues of conversely, (n-1) values of £ will produce singular- 
ities for each value of to . 

n 

Bowie illustrated this principle with two examples , which 
will be reproduced at this point: 

Example I 

Consider the characteristic equation 
f ( s ) = 0.04s 4 + 0.34s 3 + (0.2a+1.12)s 2 + (0 . 5a+B+l . 7) s 

+ 23 + 1 = 0. (2-3) 
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For ease of coefficient identification, display the terms as 




Then, if £ = 0.5, 

B. = 0 . 2u) 2 C. = -2 

In 1 

2 

B 0 = -0.5o) + 0.2oa C 0 = -a) ; 

z n n z n 



and 



B 1 C 2 B 2 C 1 “ 



-0 . 2o) 3 + 0 . 4a) 2 - to = 0, 
n n n 



Reducing , 



2 

03 - 2o3 +5 = 0; 

n n 



(2-4) 



(2-5) 

( 2 - 6 ) 



and 

03 n = 1 ± j2 . (2-7) 

Since no real values of 03 n exist, no singular line will exist 
for this characteristic equation for the value £ = 0.5. 
Example II 

f ( s ) = s 6 + 80s 5 + (20a+1600) s 4 + 840as 3 + (1600a+4006) s 2 

+ 16008s + 16008 = 0. (2-8) 
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Displaying terms , 



k 


6 5 


4 


3 


2 10 




b k 




20 : 


840 


1600 




c k 

f k 

For £ = 0 


1 80 
• 5, 


1600 


- 


400 1600 1600 




B. 


= 1600(a) 2 • 


- 840(a) 3 


C. = 


-1600 + 400a) 2 




1 


n 


n 


1 


n 


(2-9) 


B_ 


= 1600(a) 2 • 


4 

- 20(a) 


C~ = 


-1600(a) + 400a) 2 




2 

and 


n 


n 


2 


n n 




Solving , 


4 3 

go - 42go 
n n 


+ 164(a) 2 
n 


- 320o) + 320 = 0. 

n 


(2-10) 


(A) 


n = 2 ' 1115 


, 37.889 


, and 


a complex pair. 


(2-11) 



Therefore, for £ = 0.5, there exist two values of co —2.1115 
and 37.889, which will produce singular lines for equation 2-8. 

Once the conditions for singularity have been met, the 
singular line may be represented by the equation 



6 






( 2 - 12 ) 



Cadena [6] investigated the singular line when the char- 
acteristic equation coefficients are functions of three var- 
iables . 

Since the solution of equations 1-41 and 1-42 is a line 
in three dimensional space for a £-(a) pair, the definition 
singular line may be loosely applied. Where that line inter- 
sects the a-8-co-ordinate plane is the normal two dimensional 
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solution (or singular point) of equations 1-38 and 1-39, since 
y = 0. This actually reduces to his "degree of freedom" con- 
cept, whereby two roots may be specified by specifying the 
proper values of two of the variable parameters. 

A special case of the three dimensional singular line is 
when it happens to lie in one of the co-ordinate planes or in 
a plane parallel to one of the co-ordinate planes. When this 
occurs, the singular line of equation 2-12 exists for the 
specific value of y. 

An illustrative example of the foregoing was given by 
Cadena : 

Example III 
Given 

f ( s ) = s 4 + (10a+10) s 3 + (40a+56+10y+30) s 2 + ( 80a+106+252) s 

+ 70a + 106 + 25y + 340 = 0, (2-13) 

it is desired that a pair of complex conjugate roots be lo- 
cated at C = 0.5, = 2. Division of equation 2-13 by the 

2 

product of the desired roots, s + 2s + 4 , yields a reduced 
characteristic equation in a, 6, y and a remainder which must 
equal zero for the division to be exact. 

( -20y+200 ) s + (-10a -1 06-15 y+ 300) = 0 (2-14) 

or 

R l s + R 0 = 0 * (2-15) 



Obviously, , and R Q must simultaneously be zero, so 



and 



y - 10 = 0 



(2-16) 



10a + 106 + 15y - 300 = 0. 



(2-17) 
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This yields the singular line 

a + 3 - 15 = 0 (2-18) 

in the plane y = 10. 

Perhaps a more illuminating manner in which to view this 
problem is to display the coefficients in tabular form as was 
done for the two-dimensional examples. Thus 



Then, 



and 





k 


4 


3 


2 


1 


0 






b k 


- 


10 


40 


80 


70 






c k 


- 


- 


5 


10 


10 






d k 


- 


- 


10 


- 


25 






M-l 


1 


10 


30 


252 


340 




for £ = 


0.5 




= 2, 










B 1 = 


10 


C 1 


= 10 


II 

i—l 

a 


15 


II 

i— 1 

Pn 


-300 


B 2 = 


0 


C 2 


= 0 


°2 = 


40 


F 2 = 


-400 




10a + 


103 


+ 15y 


- 300 


= 0 










40y - 


400 = 


0. 







(2-19) 



( 2 - 20 ) 

( 2 - 21 ) 



This simultaneous system is the singular line formed by the 
intersection of the two planes described by equations 2-20 and 



2-21. The singular line is 

a + 3 - 15 = 0 



( 2 - 22 ) 



in the plane 



y = 10. 



(2-23) 
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Now, if = 10, 



D 2 = 20 ; 


(2-24) 


and 




10a + 106 + 15y - 300 = 0 


(2-25) 


20y - 400 = 0. 


(2-26) 


The singular line is then 




a + 6 = 0 


(2-27) 


in the plane 




• 

o 

CM 

II 

>- 


(2-28) 


If now d^ = d^ = d^ = 10, 




D x =95 D 2 = 180; 


(2-29) 


10a + 106 + 95y - 300 = 0 


(2-30) 


180y - 400 = 0; 


(2-31) 


and the singular line becomes 




a + B - — = 0 


(2-32) 


in the plane 

20 

Y = -9 • 


(2-33) 


By now it is apparent that any variation of coefficients d^ 


will result in different singular lines in 


different planes 


parallel to the a -6 plane. 




What then happens if coefficients 


and/or c^ are varied? 


Let c^ = c 4 = 10. Then 




C 1 = -70 C 2 = 160 


(2-34) 


and 




10a - 706 + 15y - 300 = 0 


(2-35) 


-1606 + 40y - 400 = 0 . 


(2-36) 
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The singular line for 5 = 0 . 5 , oo n = 2 now exists in a, 6, y- 
space at the intersection of the two planes described by equa- 
tions 2-35 and 2-36. Of what use is this information? If 
equation 2-36 is solved for y with the result then substituted 
into equation 2-35, the singular line is then "mapped" into 
the a-8 plane as 





a - 8 - 15 = 0 


(2-37) 


via 


the "transform" 






y = 46 + 10. 


(2-38) 


What 


. does this mean? For any consistent three dimensional se 


of equations 1-41 and 1-42, a two dimensional singular line 


can 


be realized by placing an appropriate 


constraint upon the 


third parameter. 






For further examination, let c 2 = 0. 


Then 




C 1 = -10 C 2 = -20; 


(2-39) 




10a - 106 + 15y - 300 = 0 


(2-40) 




- 206 + 40y - 400 = 0; 


(2-41) 


and 








4a - 6 - 60 = 0 


(2-42) 


if 








Y = j6 + 10. 


(2-43) 




Let b^ = 10. Then 






B 2 = 160; 


(2-44) 




10a + 108 + 15y - 300 = 0 


(2-45) 




160a + 40y - 400 = 0; 


(2-46) 


and 








50a - 106 + 150 = 0 


(2-47) 
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if 

y = -4a + 10. (2-48) 

What has been demonstrated in the preceding examples is 
the existence of a powerful tool for the synthesis of self- 
adaptive systems or, for that matter, any system for which a 
characteristic response is desired. How this tool is used 
will be the subject of further discussion later in this thesis. 
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CHAPTER III 



N-DIMENSIONAL PARAMETER SPACE 



Consider a set of n consistent, independent simultaneous 
equations where n = N : 

B^a + C^B + D^y t • • • + N^v + P^ = 0 

B2OI + C28 + D2Y + E26 + ... + N2V + P2 = 0 

(3-1) 



Bot + CB + DY"l'E6 + ...+Nv + P = 0. 
n n n ' n n n 

The solution of the system of equations 3-1 is the set of 

values obtained for the variables a, 6, y, 6, ..., v, which 

defines a point, which is a zero (n-n) dimensional entity, in 

n-dimensional space. 

If now only (n-1) independent equations in n = N var- 
iables can be obtained, the set is reduced to 

B^a + C^B + D^y + E^6 + . . . + N^v + P^ = 0 

B2OI + C2B + ^2^ ^* 2 ^ ■*■•••+ N2V + P2 = 0 

(3-2) 



B . a 
n-1 



+ c „-i 6 



+ D . Y 
n-1 ' 



+ E .6 
n-1 



N .v 
n-1 



+ P . = 0, 
n-1 



The solution to set 3-2 is a one [n-(n-l)] dimensional line 
defined by an arbitrary pair of the N variables. 

Similarly, for (n-2) independent equations, the solution 
becomes a two [n-(n-2)] dimensional plane defined by an 
arbitrary trio of the N variables. 
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By induction it is clear that if only two independent 
equations describe a system of N variables, there results equa- 
tions 

+ C 1 6 + DjY + E + ... + N^v + P = 0 (3-3) 

B 2 a + C 2 3 + D 2 y + E 2 6 + ... + N 2 v + P 2 = 0 (3-4) 

whose solution is an (n-2) dimensional hyperplane defined by 
an arbitrary (n-1). set of the N variables. 

If the coefficients B, C, D, E, . .., N, P are now func- 
tions of some constant pair (?,to n ), the solution to equations 
3-3 and 3-4 is an (n-2) dimensional hyperplane of constant £ , 
co n . If, when all but two arbitrary variables are set to zero 
and equations 3-3 and 3-4 remain consistent and independent, 
the point of intersection of the (a-2) dimensional constant £, 
CJ n hyperplane with the plane of the two remaining co-ordinates 
is determined. If equations 3-3 and 3-4 become dependent, the 
intersection is a line; if consistent, no intersection exists. 

For non-linear combinations of the N variables, the (n-2) 
dimensional hyperplane then becomes hyperbolic, parabolic, 
cubic, quartic, etc., in nature, which, in the independent and 
consistent case described in the preceding paragraph, gives 
rise to the possibility of multiple intersections with the two- 
dimensional co-ordinate planes. 

Although extremely little of practical engineering value 
can be extracted from the preceding discussion, the concepts 
expressed are valid and are useful for an appreciation of the 
two parameter problem. Since mental imagery extends at best 
to a relatively poor perception of three dimensions, the root 
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finding problem is limited at the present state of the art to 
the three dimensional case demonstrated by Cadena [6]. 

However, one limited use is the determination of necessary 
values of two parameters, say a and 6, given arbitrary values 
of the remaining n parameters, for obtaining a desired C , , oo n 
pair . 

Given a system described by the characteristic equation 

n , 

f ( s ) = l a, s = 0 (3-5) 

k=0 



where 

a k = b k a + c k S + d k Y + e k 5 + **• + p k' (3-6) 

one arrives at equations 3-3 and 3-4, where all capital let- 
ter coefficients are functions of the desired r , , pair. If 
variables y, 6, ..., v can be fixed or measured, the problem 

then becomes two dimensional in a and leading to solution 
by the methods outlined in chapter one. This is essentially 
what Cadena achieves by graphical means by folding the y-co- 
ordinate into the a-B plane. 

Figure 3.1 is illustrative of the three dimensional linear 
case. For the system 

B 1 a + C 1 B + D^y + F x = 0 (3-7) 

B 2 a + C 2 B + D 2 y + F 2 = 0 , (3-8) 

the constant 5, to n [{n=3}-2] hyperplane is a line in three 
dimensional space. 

For the case of non-linear combinations of variables, 
the C, lines become curves in space with the possibility 
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> 



Fig. 3.1 Intersection of Two Planes in Three Parameter 
Space 
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of multiple intersections with constant y-planes as illustrated 
in figure 3.2 for the system 

B^ + C 1 B + D 1 y + H x = 0 (3-9) 

B 2 a + C 2 B + D 2 y + E 2 a 2 + F 2 B 2 + G 2 y 2 + H 2 = 0. (3-10) 

For the condition when equations 3-3 and 3-4 are depend- 
ent, the solution is an (n-1) dimensional hyperplane described 
by either of the two "singular" equations. The same line of 
reasoning applies as in the "non-singular" case, resulting in 
the singular plane for the three parameter problem and the 
singular line for the two parameter problem, which is illus- 
trated in figure 3.3. 

Although the singular line is the most useful of all the 
"singular" hyperplanes because of the availability of graphical 
means for determining the remaining roots of a characteristic 
equation, singular theory can be readily applied to any multi- 
variate self-adaptive system or any other system where a con- 
stant complex root pair is desired. For any system with n 
parameters, variation in any or all of (n-1) parameters may be 
compensated by adjustment of the nth parameter. 
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Fig. 3.2 Non-linear Intersection in Three Parameter Space 
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Fig. 3.3 Singular Plane 
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CHAPTER IV 



INTRODUCTION TO COMPENSATION 

Since the singular line can be an extremely useful system 
characteristic, the ability to design compensation to produce 
a singular line for a specific complex conjugate root pair be- 
comes very desirable. The simplest mathematical scheme which 



comes 


to mind is that, for a system 






B^a + C 1 6 + F 1 = 0 


(4-1) 




B 2 CX + C 2 & + F 2 = 0 , 


(4-2) 


a singular line is present if , c i = C 2 ' 


and F^ = F 


Using 


the general coefficient X, if 






X 1 = X 2 , 


(4-3) 


then 


X L - x 2 = 0 , 


(4-4) 


and 


J 0 ( - 1 )kx k w >k-l (?) -v ?)] = °' 


(4-5) 



which, when expanded, becomes 

2 2 3 

- x Q + 0 ) n x 1 + (1-2?) w n x 2 + (-1-2C+4? ) w n x 3 

+ (-1+4?+4? 2 -8? 3 )o)^x 4 + (1+4?-12? 2 -8? 3 +16? 4 )w 3 x 5 

+ (1-6?-12? 2 +32? 3 +16? 4 -32? 5 )u^x 6 

+ (-1-6?+24? 2 +32? 3 -80? 4 -32? 5 +64? 6 )oo^x 7 

+ (-l+8?+24? 2 -80? 3 -80? 4 +192? 5 +64? 6 -128? 7 ) w 8 x 0 

n o 

+ (l+8?-40? 2 -80? 3 +240? 4 +192? 5 -44 8? 6 -12 8? 7 +2 56? 8 )a) 9 x o 

n y 

+ (l-10?-40? 2 +160? 3 +240? 4 -672? 5 -448? 6 +1024? 7 +256? 8 

-512? 9 ) J; 0 x 10 + . . . = 0 . (4-6) 
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If the coefficients of ? m for each term k are displayed 
as in Table 4-1, the result, neglecting the k = 0 row, which 
is (-1) for m = 0, is a lower left diagonal matrix with each 
coefficient repeating once in the same column. The repeating 
coefficient pairs alternate in sign as each column is formed. 
The values of the main diagonal-pair are (-2) m , and the suc- 
ceeding diagonal-pairs are formed by entering Table 4-2 with 
arguments 1 and m to obtain the necessary multipliers of the 
magnitude of the elements in the main diagonal-pair. 

Table 4-2 is easily formed by first making elements 1 q , 
itk and lj_ , mg equal to one ; then 



1., m. = 1., m. . + 1. ., m. 
i' j i' j-l i-l' ] 



(4-7) 



Therefore, by establishing the preceding conditions for 
all coefficients X of equations 4-1 and 4-2, compensation to 
produce a singular line can be effected. 

Because plant gain is usually the available control to 
adjust, it is convenient to equate all terms x n in the form 



Xq = + (l-2c)co n x 2 + ... , 



(4-8) 



since gain for a type-one unity-feedback plant with no zeros, 
which can always be realized mathematically by block diagram 
algebra, is the sum of all the Xq terms. 

This compensation can be easily achieved with a cascaded 
component whose transfer function numerator contains the term 
b^a + Cq 3 + f^ and whose denominator cancels the numerator of 
the uncompensated plant. 
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Table 4-1 

Table of Multipliers 



\rn 

f\ 


0 


/ 


2 


3 


4 


5 


6 


7 


8 


9 


o 


/ 


/ 


) 


i 


i 


t 


1 


1 


/ 


/ 


< 


i 


2 


3 


4 


5 


6 


7 


8 


9 


/o 


2 


t 


3 




IO 


\5 


Zt 


28 


36 


45 


36' 


3 


i 


4- 


10 


Zo 


35 


56 


84. 


/20 


IGZ 


220 


4 


_J 


—31- 


15 


36' 


TO 


tzc 


2/o 


230 


3 96* 


£ 14' 



Table 4-2 
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Example IV 



Example I in chapter II presents a characteristic equation 

0.04s 4 + 0.34s 3 + (0.2ot+1.12) s 2 + ( 0 . 5a+8+l . 7) s + 23 + 1 = 0 

(4-9) 

which has been shown to contain no singular lines. Suppose, 
for example, it is desired to have a singular line for the 
roots £ = 0.5, oj =1. If equation 4-9 represents the system 
in figure 4.1a, addition of a cascade compensator with trans- 
fer function 



r _ ba+cS+f 
G c (s) ' TE+l 

results in a new characteristic equation 

0.04s 4 + 0.34s 3 + (0.2a+1.12) s 2 + ( 0 . 5a+B+l . 7) s 

+ ba +cB + f = 0, 

where, if the system is made singular by definition, 
B 1 = -b + 0 . 2 = B 2 = -0 . 3 

C i = -c = c 2 = _1 



F 1 = -f + 0.78 = F 2 = -0.62; 



( 4 - 10 ) 



(4-lD 



(4-12) 

(4-13) 

(4-14) 



and the system then has a singular line solution for the roots 
£ = 0.5, to =1 when b = 0.5, c = 1, and f = 1.4. The com- 
pensated system is then as shown in figure 4.1b; and the sin- 
gular line for the desired roots is 

0.3a + 8 + 0.62 = 0. (4-15) 

Although, as stated in chapter two, a system may possess 
an infinite number of singular lines, one for a desired root 
pair may still not exist without compensation. 
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(b) 



Fig. 4.1 Example IV-4th Order Plant 
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Example V 



Figure 4.2a shows a system which was contrived to have 
singular lines. The characteristic equation is 



where 



s 3 + 


as + (K^ 


+8K 1 K 2 6)s + 20K 1 K 2 8 = 0, 


(4-16) 


II 

1 — 1 
PQ 


2 

0) 

n 


B ? = 

Z n 


(4-17) 


II 

i — 1 

O 


-20K 1 K 2 


C 0 = -8K, K„o) 

2 1 2 n 


(4-18) 


II 

i — 1 

Pm 


-2C«j 3 

n 


3 2 3 

F 0 = to - 45 o) - K, 0 ) . 
2 n n In 


(4-19) 



By the definition for a singular line, 



B 1 


i — 1 

u 


II 


" F 1 


1 — 1 
U 


II 


i — 1 

CQ 


- F 1 


B 2 


C 2 




1 

ho 


C 2 




B 2 


-F 2 



= 0 



Substituting , 



00 


-20K, K_ 


n 


1 2 


2 5o) 2 

n 


-8 K;l K 2 



= 4 OK, K 0 £(jj 2 - 8K.K 0 U) 2 = 0 
1 2^ n 1 2 n 



2 5u) 



n 



-2°KjK 2 



(4C 2 aj2-uT?+K .01 > -8K.K- 

n n 1 n 1 l 



(4-20) 



(4-21) 



= 80K. K 0 5^o) 3 - 20K,K„U) 
1 2 ^ n 1 ° 



3 

r'2“n 



+ 20K?K„u) - 16 K,K o 50) 3 n 

1 2 n 1 2^ n = 0 



(4-22) 



GO 



n 



2 ^ (jo. 



2 £ 0 ). 



n 



n 



, A -2 3 3 , T _ v 

(4? a) -a) +K-o) ) 
n n In 



■ 4?2 “n - “n + K l“n ‘ 4?2 “n “ °’ 



Solving, equation 4-21 gives 

C = 0.2, 

and equations 4-22 and 4-23 give 



o; = K, . 
n 1 



(4-23) 



(4-24) 



(4-25) 
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(a) 



:r\ * 


8K, Kz U+2.S) 






S(S Z 4-cK 5 -t K') 













(b) 




(c) 



Fig. 4.2 Example V-3rd Order Plant 
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The singular line for roots £ = 0.2, w n = K ]_ is described by 
both 

a - 20K 2 B - 0.4/K^ = 0 (4-26) 

from the B^ , , F^ coefficients and by 

0.4a - 8K 2 vi^B - 0.16 /k^ = 0 (4-27) 



from the B 2 , C 2 , F 2 coefficients. Dividing equation 4-27 by 
0.4, 



a - 20K 2 - 0.4 >/k^ = 0. (4-28) 

Since equations 4-26 and 4-28 must be identical, it is ob- 
vious that 

/kT = K. = co 2 = oj =1. (4-29) 

linn 

This is fortunate if the roots £ = 0.2, w = 1 are de- 
sired, but suppose a singular line is wanted for the roots 
S = 0.5, w n = 1 . 

By block diagram algebra the system can be reduced to 
that shown in figure 4.2b. Selecting a cascade compensator 
which cancels the plant numerator 



8K. K ? Bs+ba+cB + f 

G c (s) = 8K 1 K 2 B“(s+2.5) (4-30) 

where the s"*" numerator coefficient must contain B in order 
that C 2 be non-zero and where 8K 1 K 2 is inserted with the idea 
of hopefully preserving a resemblance of the original gain, a 
new characteristic equation 



s 3 + as 2 + (K 1 + 8K 1 K 2 B)s + ba + cB + f = 0 (4-31) 



is formed. For singularity 
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B 1 = 


-b+1 = 


B 2 m 1: 


b = 0 


(4-32) 


C 1 = 


-c = 


C 2 = -8 Kl K 2 ; 


o 

II 

00 

* 

H 

* 

to 


(4-33) 


F 1 = 


-f-1 = 


to 

II 

1 

H 


f = K.,-1 . 


(4-34) 



If now = K 2 = 1, the compensator becomes the readily real- 
izable 

G c (s) = s+2 .5 ' (4-35) 

with the compensated system in figure 4.2c having the singular 
line 

a - 86 - 1 = 0. (4-36 ) 

Figure 4.3 shows the parameter plane plot of the charac- 
teristic equation 

s 3 + as 2 + ( 86+1) s + 86 = 0 (4-37 ) 

with the singular line of equation 4-36 . The intersection of 
the singular line with the real root lines then determines 
the complete solution of equation 4-37 for selected singular 
line values of a and 6- It is interesting to note that the 
plots of the constant £ and constant lines form saddles as 
they approach the singular values. This is due to computer 
truncation as the program solution approaches the indeterminate 
condition . 

Example VI 

For the plant shown in figure 4.4a with characteristic 
equation 

s 3 + (a+2)s 2 + 2as + KK^ = 0, (4-38) 

cascade compensation is desired to produce a singular line 
for roots c; = 0.5, oo n = 1. The first step is to cancel the 
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Scale: a=5/in. ,0=l/in. 

Fig. 4.3 Example V-Charactenstic Equation 
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/3[s + i*L(2H,-i)«-K a -i]/3j 
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a 


KK,(s+Kz) 

















(b) 




Fig. 4.4 Example VI-3rd Order Plant 
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plant numerator while concurrently introducing the second pa- 
rameter in at least two different characteristic equation co- 
efficients so that neither nor C 2 are zero. This can be 
accomplished in at least two ways : 

Bs+ba+cB+f 



G c (s) KK. (s+K.T 



(4-39) 



or 



_ , N _ s+ba+cB+daB+f 

G C (S) kkJJ+T) ' 



(4-40 ) 



where the aB term in the numerator is needed because of the 
aB term produced by the denominator. The first method pro- 
duces the compensated characteristic equation 
4 . . . T , . „ x 3 



with 



’1 

: 1 



resulting in the singular line 



[ (K, 


,+2 ) a+2K 2 


]s 2 + 


(2K 2 a+6) s 


+ ba +c8 + 


f 










(4- 


41) 


-b + 


1 + k 2 = 


B 2 = 


2 - K 2 


(4- 


42) 


-c 


= 


C 2 = 


-1 


(4- 


43) 


-f + 


K 2 - 2 - 


II 

CN 


2K 2 - 1 , 


(4- 


44) 



(2-K 2 )a + B + 1 - 2K 2 = 0. 



(4-45) 



The compensator is then 



G c (s) - 



B { s+l+ [(2K 2 -l)a-K 2 ~l]/B} 
KK X (s+K 2 ) 



(4-46) 



which is shown in figure 4.4b. 

With the second method the characteristic equation becomes 

s^ + (a+B+2) s^ + (2a+28+a8)s^ + (2a8+l)s + ba + cB 



+ daB + f = 0, (4-47) 
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where 



making 



B 1 = 


-b + l = b 2 = 


2 


(4-48) 


o 

II 


-0+1 = C 2 = 


2 


(4-49) 


D i = 


-d+1 = D 2 = 


-1 


(4-50) 


1— 1 

II 


-f-2 = F 2 = 


- 2 , 


(4-51) 



G c ( s) 



s-a-3+2a3 
KK^Ti+FT" ' 



(4-52) 



which is shown in figure 4.4c. The singular line is 

2a + 23 - a3 - 2 = 0, (4-53) 

which may be written 

B = -^ 5 - + 2 (4-54) 

a-z 

and is plotted in figure 4.5. Since 

lim 3=2, (4-55) 

a -* 00 

3 could be set to 2 to give singular line performance for 
small variations about a large value of a. 

The preceding has demonstrated a rather simple mathemat- 
ical approach to singular line compensation, but it is suffi- 
cient both for a practical solution to the singular line 
problem and for an appreciation of a more elegant method which 
is the subject of the next chapter. 
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CHAPTER V 



COMPENSATING TECHNIQUES 



Consider the uncompensated two parameter control system 
described by 



B^a + C ^ + F^ = 0 



B 2 a + C 2 B + F 2 = 0. 



(5-1) 

(5-2) 



If this system is now compensated by some function of s and 
y, where y is a function of a and 8, equations 5-1 and 5-2 be- 
come 

B 1 a + C 1 8 + D 1 y(a,8) + F x = 0 (5-3) 

B 2 a + C 2 8 + D 2 y ( a , 8) + F 2 = 0 (5-4) 

for the compensated system. After equation 5-4 is multiplied 
by D^/D 2 , the set becomes 

B^ + C x 8 + D 1 Y(a,6) + F 1 = 0 



D 1 F 2 



D 1 B 2 , D 1 C 2 e , ^ , 01 , x A 

_ a + _ 6 + D lY (a f 8) + — = 0. 
Adding equations 5-5 and 5-6 gives 



(5-5) 

(5-6) 



D. B 



D. C 



D, F 



(B 1 + -±-p)a + ^ C i + + 2 d iY (ot / 8 ) + F ] _ + -i— = 0, 



(5-7) 



which, when solved for y(a ,&) , yields 

B, B„ C, C 



1 “l a 2 ''i '-'2 ^1 z 

Y(a ,8) = - l[ (5 i + ^)a + (^ + dJ )8 + V' (5 ' 8) 

the value of y for the compensator. Also, 
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- D 2 B i- D i B 2 



D 2 D l“ D l C 2, „ D 2 F 1 D 1 F 2 



2D, 



-)6 - 



2D 2 2D 2 ' 



D, Y (a, 6) = ( ^ ) a + ( 

x zu 2 “ LJ 2 

(5-9) 

which, when substituted into equations 5-5 and 5-6, gives 



D-B, D.. B_ D-C- D 1 C~ D-F, 

(B n - - n^-) a + (C, - 4-^ - 3 + F, - 1 



1 2D 2 2D 2 



'1 2D 2 2D 2 



1 2D, 



D 1 F 2 

2D~ 



= 0 



(5-10) 



D 1 B 2 



D 2 B 1 

2D 0 



°1 B 2 v , , °1 C 2 
2D,, ^ a + ( D„ 



°2 C 1 

2D„ 



D 1 C 2 D 1 F 2 
2D~) 3 + “dT" 



D 2 F 1 

2D„ 



D 1 F 2 

JdZ~ 



= o. 



(5-11) 



Reducing, equations 5-10 and 5-11 become 

(D 2 B 1 -D 1 B 2 )a + (D 2 C 1 -D 1 C 2 )6 + D 2 F ;l - D^ = 0 (5-12) 

(D 1 B 2 -D 2 B 1 )a + (D 1 D 2 -D 2 C 1 )6 + D.^ - D 2 F X = 0, (5-13) 



which, differing by a constant factor of (-1) , are obviously 
singular and are indeed the equation of a singular line for 
the compensated system. 

Several examples utilizing equation 5-8 to design a com- 
pensator to produce a singular line described by equations 
5-12 and 5-13 will be given for a variety of systems. 

Example VII 

Type 0-3rd Order-Feedback Compensation-One Parameter in Plant 
Denominator 

For the system shown in figure 5.1a it is desired to obtain 
feedback compensation to produce a singular line for the roots 
£ = 0.5, u> n = 1. The first step is obviously to introduce the 
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6' f 4- 

(S + 3)(S + S) (S+<x) 




(b) 




Fig. 5.1 Example VII-Type 0, 3rd Order Plant 
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second parameter. A slight amount of caution must be exercise( 
to insure that the second parameter occurs in coefficients 
other than those which will contain the compensating function 
y(a,B) so that the B coefficient in the singular line equation 
5-12 will be non-zero. This is shown in figures 5.1b and 5.1c 
which features the introduction of the compensating function. 
The resultant characteristic equation is then 

s 3 + (a+B+ 8 )s 2 + (8a+4B+y+16) s + 15a + 4y + 4 = 0 , (5-14) 

with 

B 1 = -14 C 1 = 1 D l = -4 F 1 = 3 (5-15) 

b 2 = - 7 C 2 = -3 D 2 = -1 F 2 = -8 . (5-16) 

Substitution into equations 5-12 and 5-8 gives the singular 
line 

14a + 13B + 35 = 0 , (5-17) 

produced when the system is feedback compensated by 

y (a, B) = - 4 — a - g— B - g— . (5-18) 

The characteristic equation is then 

s 3 + (a+B+ 8 )s 2 + (j^-a + ^pB + ^p) s - 6 a - p=-B - p 1 = 0 , 

(5-19) 

2 

which, when divided by s + s + 1 , the product of the constant 
roots, results in the reduced characteristic equation, or in 
this case, the remaining root, 

s+a+B+7=0. (5-20) 

Obviously, for a stable system, 

-a-B-7<0. (5-21) 
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Since, from equation 5-17, 

14 15 

8 - - n° - if - (5 - 22) 

jjo - || < 0 , (5-23) 

and 

a < 56 (5-24) 

for a stable system operated on the singular line. 

If Routh's stability criterion is applied to equation 
5-19, a stability area in the parameter plane is established 
by the boundaries formed by curves 

a + 6 + 8 = 0 (5-25) 

22a 2 + 4 3aB + 227a + 216 2 + 2236 + 708 = 0 (5-26) 

12a + 116 + 21 = 0. (5-27) 

The singular line intersects the boundary of equation 5-27 at 
a = 56 . 

The parameter plane plot of equation 5-19 is shown in 
figure 5.2, where the real root lines have a finite error in 
slope because an unresolved "bug" in the computer program for 
graphing real root lines. 

Example VIII 

Type 0-4th Order-Cascade Compensation-Two Parameters in Plant 
Denominator 

For the system in figure 5.3a, cascade compensation is 
added as shown in figure 5.3b, resulting in the characteristic 
equation 
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Fig. 5.3 Example VIII-Type 0, 4th Order Plant 
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s 5 + (a+6+k 2 k 1 +3) s 4 + [ (k 2 k 1 +3) (a+6 ) + aB + 3k^k 2 ]s 3 

+ [ (3k 2 k 1 +2) (a+B) + (^^+3) aB] s 2 + [2k 2 k 1 (a+3) (5-28) 

+ ( 3k 2 k^+2 ) aB+y ] s + 2k 2 k^aB + k^y = 0 , 
giving, for roots y = 0.5, a) n = 1, 

B 1 = c i = 2k 2 k 1 -l D = -k 2 k 1 +2 E 1 = -k ± F ± = -3k 2 k 2 +l 

(5-29) 

B 2 = C 2 = k 2 k^+l D 2 = -2k 2 k^+l E 2 = -1 F 2 = _k 2 k l _3 * 

(5-30) 

For the case of the aB product, equation 5-12 can be extended 
to 

(E 2 B 1 -E 1 B 2 ) a + (E 2 C 1 -E 1 C 2 ) B + (E 2 D 1 -E 1 D 2 ) aB + F 2^1 — F ^ F 2 ~ 



which leads to the singular line 

(l-2k 2 k 1 +k 2 k 3 +k 1 ) (a+B) + (k 2 k 1 -2-2k 2 k 3 +k 1 ) aB 



(5-31) 



+ 3k 2 k^ - 1 - k 2 k^ - 2k^ = 0. 



(5-32) 



For a linear singular line, 



^2^1 " 2 - 2k 2 k^ + k^ = 0 , 



(5-33) 



which, when solved for k^, yields 



k 2 +l 2 . / 2 

k l = 4kJ" - 4kJ v/k 2 " 14k 2 +1 



(5-34) 



Setting the radical equal to zero and solving for k 2 , 



13.93 < k 2 £ 0.07 



(5-35) 



for real values of k^ . Choosing k 2 = 13.93 in order to obtain 
the greatest possible residue for the singular roots to in- 
crease their dominance makes k^ = 0.268. The singular line 
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for the compensated system is then 

a + 6 - 1.67 = 0. (5-36) 

For the compensator, equation 5-8 can be extended to 

1 B 1 B 2 C 1 C 2 D 1 D 2 Pi F 2 

y(a,3) jl (g- + gj) a + (g- + g— ) 3 + ( + a3 + + ' 



which for = 0.268 and = 13.93 gives 

[14.47(a+B)-6.49aB-21.9] (s+0.268) 



G c (s) " 



(5-37) 



(5-38) 



s+3 . 73 

Although the pole-zero configuration of the above comoen- 
sating filter is recognized as undesirable in that it spans 

9 

more than one decade, it was chosen only to keep the singular 
line linear for illustrative purposes. A more physically 
feasible filter will result in a non-linear singular curve. 

The resultant characteristic equation then becomes 

s 5 + (a+6+6 . 73) s^ + [6 . 73 (a+6) +a6+ll . 2 ] s 3 + [1 3 . 2 (a+6) +6 . 73a6 ] s 2 

+ [21.94 (a+6)+6.71a6-21.9]s + 3.87(a+6) + 5.73aB - 5.87 = 0. 

(5-39) 

The reduced characteristic equation is then 

s 3 + (a+3+5 . 73) s 2 + [5.73(a+8)+a6+4.47]s + 6.47(a+6) 

+ 5 . 7 3a6 - 10.2 = 0. (5-40) 

Since, from equation 5-36, 

6 = -a + 1.67, (5-41) 

equation 5-40 becomes 

s 3 + 7.4s 2 + (-a 2 +1.67ot+14.13)s - 5.73a 2 + 14.13a + 0.6 = 0, 

(5-42) 
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for which a Routh ' s stability analysis yields a range of ooera- 
tion 



-0.04 < a < 2.51. 



(5-43) 



Figure 5.4 shows the parameter plane plot of the reduced 
characteristic equation. For clarity figures 5.5 and 5.6 show 
only the lines £ = 0.8 and C = 0.9, respectively. From these 
plots it is clear that the second pair of complex roots shift 
continuously for extremely small changes in a and 6. 

Figure 5.7 is the parameter plane plot of characteristic 
equation 5-39. 

Example IX 

Type l-2nd Order-Feedback Compensation -No Parameters in Plant 
For the uncompensated system of figure 5.8a, the two 
parameters must be introduced as in figure 5.8b. Tachometer 
feedback is used in order to retain the system type. Next 
the singular line compensating function is introduced with 
accelerometer feedback so that the singular line coefficients 
will be non-zero. The characteristic equation is then 



(lOy+Ds 2 + (10a+106+l) s + 10 = 0 



(5-44) 



with, for C = 0.5, to 



n 



1 




D 1 =10 F 1 = -9 



(5-45) 



b 2 = c 2 = -10 d 2 = 10 f 2 = 0 . 



(5-46) 



This produces a singular line 



10a + 106 - 9 = 0 



(5-47) 



and 




(5-48) 
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Fig. 5.4 Example VIII-Reduced Characteristic Equation 
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Fig. 5.5 Example VUI-Zeta = 0.8 
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Fig. 5.6 Example VUI-Zeta = 0.9 
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Fig. 5.7 
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Fig. 5.8 Example IX-Type 1, 2nd Order Plant 
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By substitution of the singular line relationship into equation 
5-48, 

Y = TO ' (5 - 49) 

Example X 

Type l-3rd Order-Cascade Compensation-One Parameter in Plant 
Denominator 

For the uncompensated system of figure 5.9a, the second 
parameter is introduced via tachometer feedback both to pre- 
serve the system type and to avoid a 8y product when the cas- 
cade compensator is added. This is shown in figure 5.9b. 

Then the compensator is introduced at the k = 0 and k = 1 
levels as in figure 5.9c, yielding the characteristic equation 

s^ + (a+l+k 2 k^)s^ + [ (l+k 2 k^) a+ 8 +k 2 k^] s^ + [^k^ (a+ 8 ) +4y] s 

+ 4k ] _y = 0 . (5-50) 

As usual, for t, = 0.5, to = 1, 



B 1 = k 2 k l 


0 

11 


CM 

I 

II 

1 — ! 

Q 


F 1 = -1 (5-51) 


1 1 

II 

CM 


C 2 “ k 2 k l +1 


I 

II 

CM 

Q 


F 2 = k 2 k 1 -l. (5-52) 



The singular line is then 

(-4k 2 k 1 +4k 1 )a + (-4-4k 2 k^+4k 1 ) 8 + 4 + 4k 2 k^ - 4k ] _ = 0.(5-53) 

If k 2 = 10 for a one decade filter pole-zero spread, equation 
5-53 becomes 

-36^01 + (-40k^+4k 1 -4) 8 + 40kJ - 4k x + 4 = 0. (5-54) 

The compensating function is then 

y (a , 8 ) = g^-[llk 1 a+(l+k 1 -10k^)8-l+10k^-k 1 ] , (5-55) 
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Fig. 5.9 Example X-Type 1, 3rd Order Plant 
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which gives a characteristic equation 



s 4 + (a+l+10k 1 )s 3 + [ (l+lOk-^a+B+lOkJs 2 

+ KlOk, + ii)c + (5k x + + £>6 - 2 ^- + 5 kl - |ls 

llki 1^1 2 1 2 ^1 

+ — 2 — —a + (j + - 5k£)0 - j + 5k^ - = 0 (5-56) 

and a reduced characteristic equation 

s 2 + (a+10k^)s + lOk^a + B - 1 = 0. (5-57) 



If k^ = 2, which gives good dominance characteristics to the 
singular roots, equation 5-57 becomes 

s 2 + (a+20) s + 20a +3-1=0; ' (5-58) 

and the singular line becomes 

6a + 133 - 13 = 0. (5-59) 

Substitution of equation 5-59 into equation 5-58 gives 

s 2 + (a+20) s + ^j|a = 0, (5-60) 



which may be rewritten 



, 2 54\ 

a ( s + — 



= -1 



(5-61) 



s ( s+20 ) 

for root locus analysis as in figure 5.10. Obviously the two 
non-singular roots are real, and the system is stable for all 
positive values of a. 

The characteristic equation is 



s 4 + (a+21) s 3 + (21a+B+20)s 2 + (^j-a + i|-B + ^J) s 






+ 11a 



37 37 

^ = 0 



(5-62) 
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Root Locus Plot 
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for the system compensated by 



G c (s) “ 



,11 37 q . 37, , 

(_ 8 a " I6 6 + r6 } (s+2) 

s+20 



(5-63) 



Parameter plane plots of the reduced characteristic equation 
and the characteristic equation are shown in figures 5.11 and 
5.12. 



Example XI 

Type l-4th Order-Parameter in Numerator 

When a parameter appears in the numerator of an uncom- 
pensated system, the technique used in the previous examples 
must be modified to avoid an ay or By product. In chapter 
four it was shown that the plant numerator could be cancelled 
by the compensator, but in many situations this may not be 
desirable or even possible. However, there do exist systems, 
especially electrical or electronic, whose transfer function 
factors may be physically separable, so that the numerator 
parameter may be removed as a gain or filter component as 
shown in figures 5.13a and 5.13b. In this case the compensat 
ing function may then be introduced in the numerator of a par 
allel filter having the same denominator, as shown in figure 
5.13c, where the compensated characteristic equation is 
4 5 2 

s + (a+2)s + (2a+100)s + (100a+3+y)s + 2B + 0 . ly = 0 

(5-64) 

if k = 0.1, a completely random value. Then for £ = 0.5, 




B 1 - 1 


C 1 = " 2 


D 1 = -0.1 


F 1 = 9 8 


(5-65) 


B 2 = -98 


C 2 = - 1 


°2 = - 1 


F 2 = 99. 


(5-66) 
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Fig. 5.13 Example XI -Type 1, 4th Order Plant 
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The singular line is 



10.8a - 1.96 + 88.1 = 0. 



(5-67) 



The compensating function 



y (a , 6 ) = -44a - 10.56 + 539.5 



(5-68) 



results in a characteristic equation 



s 4 + (a+2) s 3 + (2a+100)s 2 + (56a-9 . 56+539 . 5) s - 4.4a 



+ 0.956 + 53.95 = 0 (5-69) 



and a reduced characteristic equation 



2 



+ (a+1) s+a+98=0 



(5-70) 



s 



which has left-half plane roots for a > -1 and real roots for 
a > 20.8. 

A root locus plot for a > 0 is shown in figure 5.14, and 
a parameter plane plot of the characteristic equation is shown 
in figure 5.15. 

Example XII 

Type 2-6th Order-Parameter in Plant Numerator 

For the system of figure 5.16a again assume it is pos- 
sible to separate the numerator parameter as in figure 5.16b. 
Then a compensation scheme could be as depicted in figure 
5.16c, where if k^ = 10, k 2 = 100, the characteristic equation 

s 7 + 40s 6 + (a+600) s 5 + (30a+4000)s 4 + ( 300a+6+y+10000) s 3 

+ (lOOOa+206+lOlOy) s 2 + (1506+10050y) s 

+ 5006 + 50000y = 0 (5-71) 

gives, for 5 = 0.5, = 1, 

B ± = 701 C ± = -481 D x = -48991 F 1 = -9440 (5-72) 

B 2 = 971 C 2 = -130 D 2 = -9040 F 2 = -3401. (5-73) 
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Fig. 5.14 Example XI-Reduced Characteristic Equation Root Locus 
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Fig. 5.16 Example XII-Type 2, 6th Order Plant 
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The singular line is 

23.2a - 3 - 40.2 = 0. (5-74) 

y ( a , 8 ) = 0.06091a - 0.0120558 - 0.28465. (5-75) 

The resultant characteristic equation is 

s 7 + 40s 6 + ( a+ 6 0 0 ) s 5 + (30a+4000)s 4 + ( 300 . 06a+0 . 98798 

+10000) s 3 + (1061. 5a+7. 82448-287. 5)s 2 
+ (612.15a+28. 8478-2861)s + 3045.5a - 102.86 
- 14230 = 0; (5-76) 

and the reduced characteristic equation is 

s 5 + 39s 4 + (a+560) s 3 + (29a+3401)s 2 + (2 70 . 06a+0 . 9 8796 

+6039 ) s + 761.44a + 6.83656 - 9728 = 0, (5-77) 

which becomes, after application of the singular line relation- 
ship , 

s 5 + 39s 4 + (a+560 ) s 3 + (29a+3401) s 2 + ( 293 . 0 3a+5999 . 3) s 

+ 920.73a - 10000 = 0, (5-78) 

for which a root locus plot is shown in figure 5.17. Parameter 
plane plots of the reduced characteristic equation and the 
characteristic equation are shown in figures 5.18 and 5.19. 

Although the preceding examples are by no means the only 
method for effecting singular line compensation, they illus- 
trate the application of a mathematically rigorous technique 
that is both physically feasible and laboriously expedient. 
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CHAPTER VI 



RECOMMENDATIONS FOR FURTHER INVESTIGATION 

This thesis in no manner exhausts the study of parameter 
space and singularity. Although the compensation techniques 
can be extended to many types of control systems, it would be 
foolish and presumptuous • to assume universal applicability. 

Solid state operational amplifiers offer an attractive 
means for synthesizing the compensators discussed in the pre- 
vious chapters, but examination of the compensator transfer 
functions reveal that truly self-adaptive systems cannot be 
realized by these means until a voltage regulated linear 

resistance device can be achieved. 

♦ 

Digital compensating devices should be investigated for 
systems where expense is not an overriding factor. 

The possibility exists that Cadena's method of folding 
the third co-ordinate into the a-6 plane to obtain a grid of 
Y values may be repeated n-2 times for n parameter problems. 

Finally, a computer program to analytically solve the 
system 

B 1 a + C^B + Djj + + ... + N.jV + p = 0 (6-1) 

B 2 a + C 2 6 + D 2 y + E 2 6 + ... + N 2 v + P 2 = 0 (6-2) 

by iterative processes should be possible. 
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APPENDIX I 



Table of Values of 


(0.5) 


k 


U k (0.5) 
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